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Abstract. 

We study the equilibrium properties of a model for a binary mixture 
of catalytically-reactive monomers adsorbed on a two-dimensional substrate 
decorated by randomly placed catalytic bonds. The interacting A and B 
monomer species undergo continuous exchanges with particle reservoirs and react 
(A + B — » 0) as soon as a pair of unlike particles appears on sites connected 
by a catalytic bond. For the case of annealed disorder in the placement of the 
catalytic bonds this model can be mapped onto a classical spin model with spin 
values S = —1, 0, +1, with effective couplings dependent on the temperature and 
on the mean density q of catalytic bonds. This allows us to exploit the mean-field 
theory developed for the latter to determine the phase diagram as a function of q 
in the (symmetric) case in which the chemical potentials of the particle reservoirs, 
as well as the A — A and B — B interactions are equal. 



PACS numbers: 68.43.-h, 68.43.De, 64.60.Cn, 03.75.Hh 

1. Introduction 

Catalytically activated reactions (CARs) involve particles which react only in the 
presence of another agent acting as a catalyst, and remain chemically inactive 
otherwise. Usually, the catalyst is part of a solid, inert substrate placed in contact 
with fluid phases of the reactants, and the reaction takes place only between particles 
adsorbed on the substrate forming a (dilute) monolayer. These processes are 
widespread in nature and used in a variety of technological and industrial applications 

m 

The work of Ziff, Gulari, and Barshad (ZGB) 2 on the "monomer-dimer" model, 
introduced as an idealized description of the process of CO oxidation on a catalytic 
surface, as well as the subsequent studies of a simpler "monomer-monomer" reaction 
model 0, represent an important step in the understanding of CARs properties 
by revealing the emergence of an essentially collective behavior in the dynamics of 
the adsorbed monolayer. On two-dimensional (2D) substrates, first- and second- 
order non-equilibrium phase transitions involving saturated, inactive phases (substrate 
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poisoning, i.e., most of the adsorption sites are occupied by same- type particles) and 
reactive steady-states have been evidenced and studied in detail |2H3HsDB]El- Most 
of these available studies pertain to idealized homogeneous substrates. 

In contrast, the equilibrium properties of the adsorbed monolayer in the case of 
CARs are much less studied and the understanding of the equilibrium state remains 
rather limited. Moreover, actual substrates are typically disordered and generically 
the catalyst is an assembly of mobile or localized catalytic sites or islands the 
recently developed artificially designed catalysts |Zj involve inert substrates which 
are decorated by catalytic particles. Theoretical studies which have addressed the 
behavior of CARs on disordered substrates have been so far focused on the effect of 
site-dependent adsorption/desorption rates because natural catalysts are, in general, 
energetically heterogeneous [HIE]; only few studies, in particular some exactly solvable 
ID models of A + A — > reactions and a Smoluchowski-type analysis of d-dimensional 
CARs ^U], have addressed the case of spatially heterogeneous catalyst distribution. 

Recently, we have presented a simple model of a monomer-monomer A + B — > 
reaction on a 2D inhomogeneous, catalyst decorated substrate, and we have shown that 
for the case of annealed disorder in the placement of the catalytic bonds the reaction 
model under study can be mapped onto the general spin 5=1 (GS1) model with 
effective, temperature dependent couplings This allows us to exploit the large 

number of results obtained for the GS1 model ^JE| m order to elucidate, within a 
mean-field description |13| . the equilibrium properties of the monolayer binary- mixture 
of reactive monomers on a 2D substrate randomly-decorated by a catalyst. 

The organization of the paper is as follows. In Sec. we briefly present the 
model for a binary-mixture monolayer with an A + B ca ^_^ st reaction on a 2D 
inhomogeneous, catalyst decorated substrate and the mapping to a GS1 model; in 
Sec. we present the mean-field (MF) approximation. Section 0] is devoted to a 
discussion of the MF phase diagram for the particular cases of a completely catalytic 
substrate, i.e., q = 1, and of an inert substrate, i.e., q = 0, respectively. In Sec. [Sjwe 
discuss, on the basis of the results for q = 0, the MF phase diagram for general values 
of q. We conclude with a brief summary of the results in Sec. El 

2. Model of a monolayer binary mixture of reactive species on a 2D 
inhomogeneous, catalyst decorated substrate 

We consider a 2D regular lattice (coordination number z) of N adsorption sites 
(Fig. [|J, which is in contact with the mixed vapor phase of A and B particles. 
The A and B particles can adsorb onto vacant sites, and can desorb back to the 
reservoir. The system is characterized by chemical potentials [ia,b maintained at 
constant values and measured relative to the binding energy of an occupied site, so 
that fJ.A,B > corresponds to a preference for adsorption. Both A and B particles have 
hard cores prohibiting double occupancy of the adsorption sites and nearest-neighbor 
(NN) attractive A — A, B — B, and A — B interactions of strengths Ja, Jb, and Jab, 
respectively. The occupation of the i-th. site is described by a "spin" variable 



We assign, at random, to some of the lattice bonds (solid lines in Fig. ^) "catalytic" 
properties such that if an A and a B particle occupy simultaneously NN sites connected 




site i occupied by A (B), 
site i empty. 



(1) 
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Figure 1. 2D lattice of adsorption sites (small grey circles) in contact with a 
mixed vapor phase. Black (white) circles denote A (B) particles, respectively, 
solid lines denote "catalytic bonds". (R): configuration in which an annihilation 
reaction (/•%) takes place; (NR.): NN pair A-B, but no reaction because there is 
no catalytic bond between these sites. 



by such a catalytic bond, they instantaneously react and desorb, and the product (AB) 
leaves the system; A and B particles occupying NN sites not connected by a catalytic 
bond harmlessly coexist, and we assume that the reverse process of a simultaneous 
adsorption of an A and a B on a catalytic bond has an extremely low probability 
and can be neglected. The "catalytic" character of the lattice bonds is described by 
variables C<ij>> where < ij > denotes a pair of NN sites i and j, 

{1, < ij > is a catalytic bond, 

n 4-v, ■ ( 2 ) 

0, otherwise, 

and we take {C<ij>} as independent, identically distributed random variables with 
the probability distribution 

e(0 = 9*(C-i) + (i-«WC)- (3) 

Note that the probability q that a given bond is catalytic equals the mean density 
of the catalytic bonds. The two limiting cases, q = and q = 1, correspond to an 
inert substrate and to a homogeneous catalytic one, respectively. We further assume 
that the condition of instantaneous reaction A + B cat f^% st together with negligible 
simultaneous adsorption of an A and a B particle on a catalytic bond is formally 
equivalent to allowing a NN A — B repulsive interaction of strength A > 1, followed 
by the limit A — > oo, for A — B pairs connected by catalytic bonds. 

As shown in Ref. 11 , in thermal equilibrium and for situations in which the 
disorder in the placement of the catalytic bonds is annealed, i.e., the partition function, 
rather than its logarithm, is averaged over the disorder, the model under study is 
mapped exactly onto that of a GS1 model. The "effective" GS1 Hamiltonian describing 
the adsorbate at temperature T is: 

N N 

n e = -jj2 wi E - c E M + ) - H E ai + A E °i ' ( 4 ) 

<ij> <ij> <ij> i=l i=l 

where the coupling constants are given explicitely by: 

Ja + Jb - 2Jab k B T k B T 
J = -. r- hi(l - q) := J — ln(l - q) , 
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v Ja + Jb + 2Jab . k B T k B T 

K = + ln(l - q) := K Q + ln(l - q) , (5) 

g _ J A — Jb tt _ MA - PB . _ /J-A + PB 

~ 4 ' ~~ 2 ' ~~ 2 ' 

and T is the temperature. In the remaining part of this paper we focus on 
the symmetric case in which the chemical potentials of the two species are equal, 
Pa = Pb := Mi (implying H = 0), and J a = Jb '■= j, (implying C = 0). This model 
reduces to the original Blume-Emery-GrifHths (BEG) model m zer0 magnetic field 
H. 



3. Mean-field approximation of the free energy. 

The mean- field analysis follows closely the presentation in Ref.^31 an( A thus here we 
only briefly outline the main steps. The starting point is the variational principle for 
the free energy F (see, e.g., Ref.[T4"]1 

F < := Tr( P n e ) + k B TTr{phxp) (6) 

where p is any trial density matrix, i.e., Tr(p) = 1; the equality holds for p = 
exp(— (3H e )/ Tr[exp(— /3H e )], where /3 _1 = k B T and Tr denotes the sum over all spin 
configurations. 

Within the mean-field approximation the trial density p is chosen from the 

N 

subspace of products of single-site densities, i.e., p = Yi Pi\ furthermore restricting 

to the case of translationally invariant states, i.e., pi being independent of i, leads to 
trial densities of the form p = p^ . Note that this last restriction implies that within 
the present approximation the analysis cannot account for the occurrence of staggered 
states (i.e., splitting in ordered sub-lattices) and the emphasis is put on the disordered 
and ordered homogeneous states. The single site density pi minimizing the functional 
subject to the constraint Tr(pi) = 1, is 

pi=exp(-0h)/Tr[exp(-0h)], 

h = - J' Max + (A - K'Q)a\ , (7) 
where J' — zJ, K' = zK, and 

M := fa) =Tr(p lCTl ), 

Q ■= (of) = Trfaof) (8) 

are the so-called magnetization M and the quadrupolar moment Q, respectively. This 
leads to the following approximation f m f of the free energy per-site: 

f mf (M, Q) = $\jh]/N = \( J'M 2 + K'Q 2 ) - ~ In [l + 2e^ A e^ cosh(/3J'M)l .(9) 

Note that in the binary mixture language the values of the magnetization and of the 
quadrupolar moment [Eq. JSJ] represent the difference and the sum (total coverage) of 
the average densities ua and n B of A and B species, respectively: 

M = riA — n B , Q = n,A+n B . (10) 

For given values of the temperature T and of the field A = — p, [Eq.JjJJ], the 
order parameters M and Q are obtained by solving Eqs.JTJ) and JSJ. The pair (M, Q) 
characterizing the state of the system is selected from the possible solutions as the 
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one which minimizes f m t in Eq.© above. Explicitely, the equations determining M 
and Q are: 

M= 2sint ^ J/M > (11) 

exp(-/3/i - PK'Q) + 2 cosh(/3 J'M) ' v ' 

Q= 2cosh(/?J'M) 

V exp(-/3/z - (3K'Q) + 2 cosh(/3 J'M) ' 1 j 

Note that there is always a solution of these equations with M = 0, i.e., a disordered 
state (or, in the language of binary mixtures, a mixed state). 

Alternatively, one may search directly for the absolute minimum of the per-site 
free-energy function with respect to M and Q. The values M and Q at the minimum 
(minima in case of phase equilibria) will define the thermodynamically stable phase(s). 
While the first formulation is more useful for analytical work, reduced to analyzing the 
number of solutions of two coupled algebraic equations, the latter is advantageous for 
numerical calculations. In the following we shall use both of them. Before proceeding 
we note that f m f(M,Q) [Eq. ©] is an even function of M (i.e., invariant under the 
change M — ► —A/), and thus in the following we shall restrict the discussion to the 
case M > 0; the states with M < are immediately obtained via a change of sign. 
This is a consequence of the symmetry in the chemical potentials (p. a — Ms) or 7 m 
the magnetic language, of a vanishing magnetic field H = 0. In other words, in the 
space spanned by (T, H) the phase diagrams in the plane H = along the H = + 
side and any equilibrium state characterized by (M > 0, Q) have corresponding phase 
diagrams and states {—M, Q) located on the H = 0~ side. 

4. Homogeneous catalytic or catalytically inert substrates. 

4-1- The case of a homogeneous catalytic substrate: q = 1. 

The case of a homogeneous, completely catalytic substrate can be studied analytically 
because of the particular form of the interaction parameters K' and J'. In the limit 
q — > 1, Eq. implies: 

T t ksT 1 \ l 

J = ZJq — Z 2 U1 (l ~ 1) K +°° 7 

K' = zJ + z^ln(l- 9)^-00, (13) 

while J' + K' = zj remains finite in this limit. The analysis of Eqs. ljTJl and (|12H 
proceeds as follows. From Eq. l|12|l . the solutions with M = 0, i.e., disorder states, 
have the quadrupolar moment 

2 

® ~ a- 1 ™ 00 cxp(-Pn) exp(-f3K'Q) + 2 ~ ° ' ^ 

for any finite temperature and finite chemical potential /x. This corresponds to an 
empty lattice state. 

Taking the ratio of Eqs. ljTJl and l|12|l we find that the ordered states, M 5^ 0, 
satisfy 

Q= lim M coth(^ J'M) = \M\, (15) 

J'— ^ + 00 
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In this limit the substrate is occupied by a single species, either A or B with equal 
probability. With K'+J' = jz, for any finite temperature and finite chemical potential 
M is determined from 

f - exp(-2/3J'M) 



M = 

exp(/3A) exp(-0jzM) + 1 + exp(-20J'M) 

J^ + oo I (lg) 



exp(— /3^) exp(— (3jzM) + 1 ' 

It is easy to see that Ea. (|l(j|) has a solution < \M\ < 1 for any T > and any /x. 
Thus, in virtue of Ea. l(T5|) . the lattice is occupied, with equal probability, either by 
\M\ x N particles of species A or by \M\ x N particles of species B. 

4-2. The case of an inert substrate: q = 0. 

In the case of a catalytically inert substrate, 9 = 0, we have 
ji j-2J AB , j + 2J AB 

J = z - = zJ , K = z = zK , (17) 

and thus the model reduces to the classical BEG model, whose mean-field 
approximation has been analyzed in detail in Ref. [Ej- I n the following we briefly 
summarize the main aspects of the phase diagram of the BEG model such that in the 
general case g^Owe can isolate the effects solely due to the disordered distribution of 
the catalyst. Moreover, as will be shown in the next section, the equilibrium properties 
of the adsorbate in the general case q ^ can be easily rationalized from the ones on 
the inert substrate. 

First we note that for non-interacting particles, i.e., for j = Jab = so that 
J' = K' = 0, Eq.jnj implies that M = 0, and thus Ea.lJT2|) leads to 

2 ,n A = n B =Q/2 = —l—, (18) 



** ~ exp(-/3 M ) + 2 1 + 2 f 

where f := exjp(Pfi) is the fugacity, rendering the classical Langmuir adsorption result. 

For the case of non-zero interaction parameters j and Jab, only the qualitative 
features of the phase-diagram can be analytically derived (for details see Ref.fT^]'): here 
we shall present phase diagrams obtained via direct numerical minimization of the free 
energy function f m f [Eq.©]. Using zJ as energy scale, the system is characterized 
by the parameter k := K'/J' — K /J , the scaled temperature (thermal energy) t, 
and the scaled chemical potential u: 

t := k B T/{zJ ) , u := ^/{zJ ) . (19) 

We shall discuss below the phase diagrams, as well as the behavior of the order 
parameters M and Q, in the u — t plane at given values of ko > 0. The reason 
for the latter restriction is that for sufficiently negative values kq < < 0, where 

(tr) 

Kq is a threshold value, it is known that the system will split into ordered sub- 
lattices ^Hl- However, such states are not captured by the present formulation of the 
mean-field equations which assume translational invariance. 

Before proceeding with the numerical analysis, we list several general features of 
the phase diagram that can be obtained from the analysis of Eas. (lll|) and p2|l. 
(i) For M/0, Ea. fTT|) can be written as 

2[tecosh(x) — sinh(x)] + tx exp(— u/t) exp(— KoQ/t) — 0, (20) 

V v ' " v ' 

:— gi(x-.t) :—g2{x;t,u) 
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where < x := M/t < 1/t. Since g\{x;t > 1) is a strictly increasing function 
of x, it follows that gi(x > 0;t > 1) > gi(Q;t > 1) = 0; moreover, one has 
g2{x > 0;t > 0,u) > 0, from which one can conclude that for t > 1 Eq.J^OJ) has 
only the solution x = 0, i.e., there are no ordered states for t > 1. 

(ii) For /i < fJ. c (T) < 0, the positive term g2(x;t,u < 0) dominates over the term 
<7i(x;i), which is bounded from below, and thus in this range Eq. 1201) has only the 
solution x — 0, i.e., there are no ordered states for u < u c (t) < 0. This is an intuitive 
result: \i < corresponds to a preference for desorption; thus at low negative chemical 
potential the substrate is covered by a low density two-dimensional gas (which is a 
mixed phase). 

(iii) For /j — > oo and finite temperatures, g2(x;t,u — » oo) — > 0, and thus Eq. i|20[l 
reduces to gi(x;t) = which is equivalent to tanh(x) = tx. It follows that for t < 1 
there is always a unique solution > 0, and therefore there is an ordered state 
M(t < l,n — > oo) 7^ such that M(t /* — > oo) — » 0, i.e., if there is a phase 
transition at t — 1 (for — > oo), it is a continuous order-disorder transition. 

(iv) In the case of disordered states (i.e., M — 0) Eq.JT^J reduces to 

2ty + tyexp(-uft) exp(-K y) -2 = 0, (21) 
■■=93(y\u,t) 

where < y := Q/t < 1/t. For /i > and finite temperatures, gz{y;u > 0,t) is a 
strictly increasing function of y, and thus Eq. I|21|l has an unique solution, i.e., in the 
region /i > there are no phase transitions between disordered states. 

We now turn to a detailed discussion of the phase diagrams. The most 
complicated behavior occurs for intermediate values of ko, e.g., no ~ 3 13 . For 
kq = 3.0, in Fig. [2] we show (color coded) the order parameters M = \tia — ng\ 
and Q = ua + ns, as well as the lines corresponding to the various phase transitions 
composing the phase diagram. 

At low temperatures, i.e. t < tp 2 ~ 0.72 and low negative chemical potential, 
i.e., u << up 2 ~ —2.0, the substrate is covered by a very low density, two-dimensional 
mixed gas (M = 0, Q -C 1), in agreement with (ii) above. Increasing the chemical 
potential u at fixed temperature t < tp 2 , the system undergoes a first-order phase 
transition at u = u c (t) ~ —2.0 (the white dashed line located at u c (t) ~ —2.0) upon 
which the density of the monolayer increases abruptly to almost one [as indicated 
by the dark blue color in Fig. |2Jb)] and the monolayer also (partially) demixes, i.e. 
< M < Q [as indicated by the lighter blue or green color in Fig. [Ha)]. Thus in 
the region u > u c (t) as anticipated in (iv) above the substrate is covered by a dense 
A-rich monolayer (for H = [ia — t l B — * + ; with equal probability, a dense £?-rich 
monolayer forms for H = \ia — Hb — > 0~). Therefore this first-order transition line, 
which is located at almost constant /i ~ —2.0 and extends from P2 to t — > 0, is a line 
of triple points. 

We focus now on the region u > u c (t). At constant u, upon approaching from 
below the line P 2 -Pi, which continues to 11 — > 00, the demixing is less and less 
pronounced. Increasing the temperature at fixed u > up 1 , upon crossing the line 
starting at Pi (solid black line in Fig. |2J the dense A-rich (or B-rich) monolayer 
undergoes a second-order phase transition (both M and Q are changing continuously 
there; note in Fig. |2{a) the thin band of yellow color, which ends at P%, corresponding 
to very small but non-zero values of M) such that at high temperatures the substrate 
is covered by a mixed dense monolayer. As discussed in (i) and (iii) above, this line 
of critical points stays below t = 1 for all values of u, and it approaches t = 1 for 
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Figure 2. (a) Magnetization M and (b) quadrupolar moment Q for an inert 
substrate (q = 0) as functions of t := ksT/ (zJo) and u := /i/(zJo) for reo = 3.0. 
The color coding (shown at the right) is the same for both figures, and linearly 
interpolates between zero (red) and one (dark blue) via yellow, green, and light 
blue. Pi is a tricritical point, P2 is a quadruple point, and P3 is a critical end 
point. Solid black lines are lines of critical points and indicate second-order phase 
transitions, while the dashed white lines are lines of triple points indicating first- 
order transitions. 



u — > 00. Upon crossing the line segment P2-P1 (white dashed line) the transition from 
the dense A-rich (or B-rich) monolayer to the dense mixed monolayer is of first order 
with a jump in both coverage and composition, i.e., in both Q and M (using the color 
code, for M this is indicated by the transition from red to green, without yellow in 
between, while for Q by the direct transition from dark blue to light-blue and green). 
The line segment P2-P1 is also a line of triple points since there three phases coexist: 
dense A-rich, dense _B-rich, and dense mixed monolayer, respectively. 

For u < u c (t), upon crossing (from below) the line P2 — P3 the system undergoes 
a first-order transition from a low density mixed monolayer to a higher density mixed 
monolayer: there is a jump in the coverage, i.e., Q varies discontinuously [as indicated 
in Fig. |U^b) by the color change from light green to light blue (without dark green 
in between)], while M remains zero. The jump in Q upon crossing the line P2 — P3 
decreases as the crossing point approaches P3, and it becomes zero at P3. Thus P3 is 
a critical point. Note that for temperatures t such that tp 2 < t < tp 3 , e.g., t = 0.73, 
increasing u at constant temperature from small negative values toward positive values 
will drive the state of the monolayer from the mixed gas phase towards the A-rich (or 
i?-rich) dense phase via two consecutive first-order phase transitions, corresponding 
to crossing the line P2 — P3 (with a jump only in Q) followed by crossing the line 
P2 — Pi (with a jump in both Q and M). 

The point Pi is a tricritical point (it belongs also to the critical lines of the A-rich 
dense phase — > mixed gas and i?-rich dense phase — > mixed gas transitions). P2 is 
a quadruple point, i.e., at P2 four phases coexist: dense A-rich, dense B-rich, dense 
mixed, and dilute mixed monolayer, respectively 

Varying k will lead to topological changes only near the points Pi, P2, P3 as 
follows. For Ko C 1, the line P2-P3 and thus the point P2 do not occur [there is 
no jump in the total coverage if the density is increased while keeping the monolayer 
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mixed, i.e., in that region which is red in Fig. Ufa)], so that at the tricritical point Pi 
the line of triple points connects directly with the one corresponding to the second- 
order phase transitions. With increasing Kq, the line P2-P3 emerges, and for k < 3 
the behavior is the same as the one at kq = 3. For kq > 3 the only change is that the 
point P3 is located at higher values of t than Pi . With increasing kq > 3 Pi is shifting 
towards Pi and eventually reaches the first-order triple line such that for large kq, 
e.g., no ^ 5, Pi (which at low and medium values of Ko is a tricritical point) merges 
with P2 forming a critical end point. 



5. The case of a disordered substrate: < q < 1. 

In the case of a disordered substrate, < q < 1, one has 
J' = zJ -^Hl-Q), 

K > = zK + ln(l - q) . (22) 

Using again zJq as the energy scale, the system is now characterized by the two 
variables t and u defined in Subsec. 14.21 a disorder parameter q defined as 

q = 1 - (1 - qf , i.e., z ln(l -q)= ln(l - q) , (23) 

and the parameter 

Here we indicated explicitely that now the ratio k depends both on the temperature 
and on the disorder parameter q. [Note that n(t, 0) = k(0, q) — kq.] Therefore, we shall 
discuss the phase diagrams, as well as the behavior of the order parameters M and 
Q, in the u — t plane for given values of kq and q. 

Ea. i|24|) implies that for any given < q < 1 and kq > 0, the ratio K(t, q) becomes 
negative at high enough temperatures, i.e., n{t > t tr ,q) < 0, where the threshold 
temperature U r is given by 

t tr {q;K ) = 2«o/mi-g)|; (25) 

note that for a given kq, i.e., for a given mixture, t tr is a decreasing function of q. As 
already mentioned, for sufficiently negative values of n the system splits into ordered 
sub-lattices [T^l, which generally leads to a significant decrease in the yield of the 
catalytic reaction. Using t tr as a measure of this tendency, Ea. (|25|l implies that it 
is desirable to run the reaction at low enough temperatures in order to maintain a 
mixed monolayer, and that this range of temperatures decreases with an increasing 
mean density of catalytic bonds. 

At constant temperature, and for fixed values of q and of kq, the parameter n(t, q) 
is constant; for two temperatures ti and £2 > ti, it satisfies n{t2 > ti,q) < n(ti,q) 
(in particular, one has n(t,q) < kq). Thus, for any chosen temperature i < t tr (q;Ko) 
the isotherms M(t,u) and Q(i,u) can be simply read as the ones corresponding to 
an inert substrate, as in Subsec. 14.21 at the same temperature t but for a binary 
mixture with kq ~ n(t,q) < kq. Since larger values of i correspond to smaller 
values of Rq, the phase diagram is expected to be more similar to one at low kq 
on an inert substrate, and thus only the first-order mixed gas — > A-rich (i?-rich) 
dense monolayer transition and the second-order transition lines joining at a tricritical 
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point are generally present. The tricritical point Pj shifts towards increasing values 
of u with increasing q, and the first-order transition line is a smooth curve (rather 
than consisting of two segments P2 — Pi and Pi — > zero temperature), which runs 
from u ~ —2.0 at low temperature (corresponding to the location in the case kq) to 
Pi. These qualitative features can be seen in Figs. [3] and 0] where we show results 
corresponding to kq = 3.0 for q = 0.3 and q — 0.6, respectively, which allows one a 
direct comparison with the phase diagram on the inert substrate. Note that in these 
cases the threshold temperatures [Ea. (|25|l ] for splitting into ordered sub-lattices are 
very high [i tr (0.3; 0.3) ~ 16.8, £ tr (0.6; 0.3) ~ 6.54], thus the present version of the 
mean-field analysis is justified in the range t < 1 we are interested in. 



(a) [M] (b) [Q] 
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U U 

Figure 3. (a) Magnetization M and (b) quadrupolar moment Q as functions 
of t := ksT I (zJo) and u := /i/(zJo) for kq = 3.0 and q = 0.3. The color coding 
(shown at the right) is the same for both figures, and linearly interpolates between 
zero (red) and one (dark blue) via yellow, green, and light blue. Solid black lines 
are lines of critical points and indicate second-order phase transitions, while the 
dashed white lines are lines of triple points indicating first-order transitions. 
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Figure 4. Same as Fig. El for q = 0.6. 
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For an efficient catalytic reaction, the system would have to be operated such 
that the substrate is covered by a mixed, not too dilute monolayer, i.e., not too low 
positive chemical potential and temperatures above the second-order transition line, 
i.e., t ~ 1.0, but below the threshold temperature tt r (q~', k-o) for which the splitting 
into ordered sublattices may occur. For a given binary mixture, i.e., given kq, the 
requirement tt r < 1 implies an upper bound q op on the mean density of catalytic 
bonds, i.e., 

t tr < 1 =>■ q < q op — 1 - exp(-2«o) . (26) 

This is somewhat surprising: a substrate which is only partially decorated by the 
catalyst, i.e., q < q op < 1, would be the optimal choice because it avoids transitions 
into a passive state (poisoned substrate). For binary mixtures with small values of Kq 
the upper bound given above implies a drastic constraint (e.g., at k = 0.5, q op ~ 0.63), 
while for binary mixtures with large values of kq, the above constraint is basically 
irrelevant (e.g., at Ko = 3.0, q op ~ 0.998). 

6. Summary 

Within a mean- field approach we have studied the equilibrium properties of a model for 
a binary mixture of catalytically-reactive monomers adsorbed on a two-dimensional 
substrate, decorated by randomly placed catalytic bonds of mean density q. Our 
analysis here has been focused on annealed disorder in the placement of the catalytic 
bonds and on the symmetric case in which the chemical potentials fiA arL d fj-B , as well 
as the interactions J a and Jb of the two species, are equal. 

We have shown that in the general case < q < 1 the mean-field phase diagram 
and the behavior of the composition ua — ub and of the total coverage ua + ns can 
be extracted from the corresponding results on an inert substrate (q — 0). We have 
determined certain restrictions on the temperature, at which the system is operated, 
as well as a somewhat surprising upper bound q op of the density of catalytic bonds, 
which have to be obeyed in order to maintain the monolayer in a mixed and thus 
active state. Even in this highly symmetric case studied here, the phase diagrams 
are rather rich, containing, e.g., lines of second-order phase transitions and tricritical 
points. We have pointed out the likely occurrence of staggered phases at intermediate 
temperatures. 

Finally, we note that in the generic case in which [iA Hb and J a ^= Jb, the 
behavior is expected to be even richer (see, e.g., Ref. (TBI for a detailed discussion of 
the phase diagrams for the general spin S = 1 model). These aspects, as well as the 
issue of staggered phases or that of quenched disorder, are left for future work. 
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